Normal forms, stability and splitting of invariant 
manifolds I. Gevrey Hamiltonians 



Abed Bounemoura * 
December 7, 2012 

Abstract 

In this paper, we give a new construction of resonant normal forms with a smaU re- 
mainder for near-integrable Hamihonians at a quasi-periodic frequency. The construction 
is based on the special case of a periodic frequency, a Diophantine result concerning the 
approximation of a vector by independent periodic vectors and a technique of composition 
of periodic averaging. It enables us to deal with non-analytic Hamiltonians, and in this 
first part we will focus on Gevrey Hamiltonians and derive normal forms with an exponen- 
tially small remainder. This extends a result which was known for analytic Hamiltonians, 
and only in the periodic case for Gevrey Hamiltonians. As applications, we obtain an 
exponentially large upper bound on the stability time for the evolution of the action vari- 
ables and an exponentially small upper bound on the splitting of invariant manifolds for 
hyperbolic tori, generalizing corresponding results for analytic Hamiltonians. 

1 Introduction and main results 

1. Let re > 2 be an integer, T" = M^/Z" and Br be the closed ball in M", centered at 
the origin, of radius R > with respect to the supremum norm. For e > 0, we consider a 
near-integrable Hamiltonian of the form 

(Hie,i) = hii) + fie,i) 

\\f\<e«l 

where (0,1) G P/j = T" x are angle- action coordinates for the integrable part h and / 
is a small perturbation in some suitable topology defined by a norm | . |. The phase space 
Vr is equipped with the symplectic structure induced by the canonical symplectic structure 
on T" X M" = T*T". For e = 0, the system is trivially integrable: the action variables 
are integrals of motion and their level sets {/ = Iq}, Iq G Bji, are invariant embedded tori 
on which the flow is linear with frequency V/i(/o). An important subject of study is the 
dynamical properties of these systems when e > 0. 

Here we will be interested in the local dynamics of the perturbed system around an 
invariant torus of the integrable system, which, without loss of generality, we assume is 
located at / = 0. The qualitative and quantitative properties of this invariant torus with a 
linear flow is then determined by its frequency vector w = V/i(0) G M". Let us say that a 
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vector subspace of M" is rational if it has a basis of vectors with rational (or equivalently, 
integer) components, and we let F = F^^ he the smallest rational subspace of M" containing 
ui. li F = M", the vector u is said to be non-resonant and the dynamics on the invariant 
torus To = {-/^ = 0} is then minimal and uniquely ergodic. If is a proper subspace of M" of 
dimension d, the vector oj is said to be resonant and d (respectively m = n — d) is the number 
of effective frequencies (respectively the multiplicity of the resonance): the invariant torus 
To is then foliated into invariant d-dimensional tori (which are just images of translates of F 
under the canonical projection vro : — t- To) on which the dynamics is again minimal and 
uniquely ergodic. The example of a resonant vector to keep in mind, and to which the general 
case can be actually reduced, is a; = {w, 0) G = M*^ x with w ^ M.'^ non-resonant, 
for which F = M.'^ x {0}. From now on, we will assume that 1 < d < n (the case d = 
corresponds to an invariant torus which consists uniquely of equilibrium solutions), and the 
particular case d = 1 will play a special role: in this case, writing uj = v, the torus is foliated 
by periodic orbits and if T denotes the minimal common period, then Tv E Z", and such a 
vector V will be called periodic. 

In order to study the local dynamics of the perturbed system around such an invariant 
torus, it is important to quantify the minimal character of this linear flow with frequency to 
(on the whole invariant torus if d = n or on each leaf if d < n). In general, to such a vector 
UJ one can associate a constant Qo; > and a real-valued function defined for all real 
numbers Q > Qi^, which is non-decreasing and unbounded, by 

"^UQ) = inax{\k-uo\-'^ I G Z" nF, < |fc| < Q} (1) 

where • denotes the Euclidean scalar product and | . | is the supremum norm for vectors (see 
for instance |BF12] for much more detailed information, where this function is denoted by 
^'^). Note that by definition, we have 

|A;-c^|>^^, keZ'^nF, o<\k\<Q. 

Let us say that a vector uj is Diophantine if there exist constants 7 > and t > d — 1 such 
that ^uj{Q) < 7~^Q^, and let us denote by r2d(7, r) the set of such vectors. For d = 1, u = v 
is T-periodic and F = F^ is just the real line generated by so that for any non-zero vector 
k G "Z" D F, as /c is collinear to v and k and Tv are non-zero integer vectors, we have 

\k.v\-^ = \k\-'^\v\-'^ = T\k\-^\Tv\-'^ < T 

hence ^v{Q) < T and therefore any T-periodic vector belongs to Oi(T~^,0). 

2. For an analytic Hamiltonian system, it is well-known that in the neighbourhood of 
an unperturbed invariant torus with a Diophantine frequency vector, the system can be 
analytically conjugated to a simpler system where the perturbation has been split into two 
parts: a resonant part, which captures the important features of the system and whose size 
is still comparable to e, and a non-resonant part, whose size can be made exponentially small 
with respect to e~", where the exponent a > depends only on the Diophantine exponent r. 
The result can also be extended to an arbitrary vector uj £ M", in which case the non-resonant 
part is exponentially small with respect to some function of e^^, this function depending only 
on (see |Boul2b ] for d = n). Such simpler systems are usually called resonant formal forms, 
and were first obtained in this context by Nekhoroshev (see |Nek77j . or |Pos93] for a nicer 
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exposition) with the aim of estabhshing exponentiahy long (and global) stability estimates for 
the evolution of the action variables, under some geometrical assumptions on the integrable 
part. Apart from deriving such stability estimates, these resonant normal forms can also be 
used to prove the existence of invariant hyperbolic objects such as tori or cylinders of lower 
dimension. For hyperbolic tori (also called whiskered tori), which by definition possess stable 
and unstable manifolds (also known as whiskers), these normal forms can be used to estimate 
the "angle" between these invariant manifolds at an intersection point: this is usually called 
the splitting of invariant manifolds (or "separatrices" ) . 

The aim of this paper is to construct resonant normal forms with a small remainder in the 
broader class of Gevrey Hamiltonians, for an arbitrary frequency vector uj. The case d = 1 is 
known and our main theorem extends this result for any \ < d < n. Our proof uses a method 
of periodic approximations: it is based on the simple case d = 1, a Diophantine result on the 
approximation of a vector by independent periodic vectors proved in |BF12j and a technique 
of composition of periodic averaging first used in |BN12j . Note that our main result answers 
a question which was asked in jLMS03j . concerning the "interaction of Gevrey conditions 
with arithmetic properties in normal forms". In a second paper |Boul2a) . we will prove the 
corresponding result for finitely differentiable Hamiltonians using a similar method (but with 
fairly different technical details). 

We give two applications of our main result. As a first application, we prove local and 
exponentially long stability estimates for the evolution of the action variables for perturbations 
of non-linear integrable systems, without any condition on the integrable part. We also prove 
global and exponentially long stability estimates for perturbations of integrable linear systems, 
and show that the latter result is essentially optimal. If one is interested in global stability 
estimates for perturbations of non-linear integrable systems, one needs to impose a geometric 
condition on the integrable part. It is known that for convex integrable systems, an essentially 
optimal result can be obtained using only the case d = 1 (see |Loc92j for analytic systems, 
|MS02] for Gevrey systems and [BoulO] for finitely differentiable systems). However, it seems 
that the general case 1 < d < n is needed if one is interested in such global estimates for 
non-convex integrable Hamiltonians. In [BN12] and [Boullj . special normal forms adapted to 
the problem were constructed using only periodic approximations: however, the quantitative 
results obtained there were very far from being optimal, and in particular our main theorem 
in this paper gives much more general and precise normal forms that should be useful in 
trying to also obtain essentially optimal results for perturbations of non-convex integrable 
systems. As a second application, we will prove an exponentially small upper bound for the 
splitting of invariant manifolds for a hyperbolic torus. 

3. Let us now state precisely our results, starting with the regularity assumption. Recall 
that n > 2 and i? > have been fixed, and that our phase space is P/j = T*^ x B^^ where 
Bn is the closed ball in M", centered at the origin, of radius R with respect to the supremum 
norm. 

Given a > 1 and L > 0, a real-valued function / G C°°{'Dr) is (a, -L)-Gevrey if, using the 
standard multi-index notation, we have 

Iflc^^H-Dn) = \f\-,L,l,R < OO, = ^'""(Z!)-"|9'/bo(l.«) (2) 

where | . Icocdu) is the usual supremum norm for functions on Dpt- The space of such functions, 
with the above norm, is a Banach space that we denote by G'^'^{V>r). In the sequel, we shall 
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forget the dependence on the domain in the notation and simply write | . \ a,L = | • \g'=''^{Vr) 
and I . \a,L,i = I • \a,L,i,R- Next, for p > 1, a vector-valued function F G C°°{Vii,W) is {a,L)- 



Gevrey if F = (/i, . . . , /p) with fi £ for 1 < i < p. The space of such functions, 

that we denote by M^), is also a Banach space with the norm 

\^\G°'''-{Vit,Rp) = l/^lG°-^'(I'J^)• 
Similarly, we shall simply write | . \a,L = I • lG"'^(X'fl,Rp)- identifying functions with values 
in T with 1-periodic real-valued functions, we can also define G'^'^{'Dji, T" x M") and | . \a,L = 
I . Ig"'^(X'h,T"xR")- Finally, we will also use the notation G° = Ul>o to denote the space 
of functions which are (a, L)-Gevrey for some L > 0, on appropriate domains. 

Note that analytic functions are a particular case of Gevrey functions, as one can check 
that G^'^{Vji) is exactly the space of bounded real-analytic functions on Vji which extend as 
bounded holomorphic functions on the complex domain 

VUVn) = {(9,1) G (C7Z") X C" I \Iie)\ < L, diI,Bn) < L}, 

where I{0) is the imaginary part of ^, | . | the supremum norm on C" and d the associated 
distance on C". 

Since we will be only interested in local properties, it will be enough to consider first linear 
integrable Hamiltonians, that is = = uj-I. As we will see, our result for an arbitrary 
non-linear integrable Hamiltonian h will be obtained from the linear case by a straightforward 
localization procedure. Therefore we shall first consider a Hamiltonian H G G'^'^{T)ji) of the 
form 

iH{e,I) = L{I) + f{0,I), (9,1) eVn, 

\\f\a,L<e. ^ 

We denote by { . , . } the Poisson bracket associated to the symplectic structure on Vji. For 
any vector w G M", let X^^ be the Hamiltonian flow of the linear integrable Hamiltonian 
lw{I) = w ■ I, and given any g G C^(Pij), we define the average (along the linear flow of 
frequency lo) of g by 

[gU = lim - /" goX^dt. (3) 

Note that {g, 1^} = if and only if goX^ = 5 if and only if 5 = [g]w Recall that the function 
has been defined in ([1]), then we define the functions 

A^(Q) = Q^UQ), Q > Qu., Kix) = sup{Q > I A^(Q) <x}, x> A^(Q^). 
Our first result is the following. 



Theorem 1.1. Let H he as in (|Gip . There exist positive constants c, ci, C2, C3, C4, C5 and 
C that depend only on n,R,uj,a and L such that if 

AUce-')>cu (4) 

then there exists a symplectic map ^ G G"'^(P/j/2) ^-r); where L = CL, such that 

Ho<^ = l^ + [f]^ + g + f, {g,L}=0 
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with the estimates 

|$-Id|„^^<C2A:(ce-i)-i (5) 

and 

\9\^^^ < C3eAUce-')-\ I/L,l < C4eexp (^-c,AUce-')^) . (6) 

On the domain the above theorem states the existence of a symplectic Gevrey 

conjugacy, close to identity, between the original Hamiltonian and a Hamiltonian which is 
the sum of the integrable part, the average of the perturbation whose size is of order e, a 
resonant part which by definition Poisson commutes with the integrable part and whose size 
is of order e(A* (ce~^))~^, and a general part whose size is now exponentially small with 
respect to A* (ce"^)^ . The first terms of this Hamiltonian, namely l^j + [f]uj + g, is what is 
called a resonant normal form, and the last term / is a "small" remainder. 

Concerning the size of this remainder, in |Boul2b] an example in the analytic case for 
d = n is given to show that this result is "essentially" optimal: we will show below in ^ 
that this example can be easily adapted to the Gevrey case and for any 1 < d < n, and that 
our estimate for the remainder in the Gevrey case is also "essentially" optimal. The word 
"essentially" should be understood as follows: given uj, one can always construct a sequence 
of positive real numbers ej = going to zero as j goes to infinity, and a sequence of 

ej-perturbations fj, such that the estimate for the remainder fj cannot be improved (see the 
next section for a precise statement). This, of course, does not preclude this estimate to be 
improved for other values of s. 

Now in the Diophantine case, the estimates of Theorem 11.11 can be made more explicit. 

Indeed, in this case we have the upper bound "^^{Q) < 7~^Q^ which gives the lower bound 

1 

A*(ce~ ) > (c7e~ ) 1+^ . The following corollary is then straightforward. 

Corollary 1.2. Let H be as in ()Gip . and uj G 0^(7, r). There exist positive constants c, ci, 
C2; C3, C4, C5 and C that depend only on n,R,uj,a and L such that if 

e < cc-^^ 7, 

then there exists a symplectic map $ G G°'^^{Vjii2,Vr), where L = CL, such that 

Ho^ = l^ + [f]^ + g + f, {g,l^} = Q 

with the estimates 

and ^ 

\g\^X ^ C3e(c"V"^e)^, l/L,L ^ C4eexp (^-C5(c7e"^) "(1+-) ) . 

Even though the size of the remainder is essentially optimal as we already explained, in the 
Diophantine case we believe that the other estimates, concerning the size of the transformation 
<I> and the resonant term can be improved. As a matter of fact, these improvements are 
known in the analytic case for d = n (see the comments at the end of this section). 

4. We can now state a local result for a perturbation of a general non-linear integrable 
Hamiltonian system, that is for a Hamiltonian H G G°''^{T>ji) of the form 

H{ej) = h{I) + f{ej), (9,1) €Vn, 

Vh{0)=UJ, \h\a,L<l, \f\a,L<e. 
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For a "small" parameter r > 0, we will focus on the domain Vj- = T'^ x B^, which is a 
neighbourhood of size r of the unperturbed torus To = x {0}. 

Since we are interested in r-dependent domains in the space of action, the estimates for the 
derivatives with respect to the actions will have different size than the one for the derivatives 
with respect to the angles. To distinguish between them, we will split multi-integers I G N^*^ 
as / = iliM) G N" X N" so that 9' = dl^d\^ and \l\ = \h\ + \l2V Let us denote by Id/ and Idg 
the identity map in respectively the action and angle space, and for a function F with values 
in T" X M", we shall write F = {Fg, F/). 

Theorem 1.3. Let H be as in ()G2p . There exist positive constants c, cq, cj, cg, cg and cio 
that depend only on n,R,uj,a and L, such that if 

V^<r<R, A*Jcr-^)>ce, (7) 
then there exists a symplectic map <I> € G°^(T>j./2:'^r) such that 

Ho<!> = h + g+[f]^ + f, {g,L} = 0. 
Moreover, for L = CL and any I = (lijh) £ N^", we have the estimates 

- Id7L,z,z < cjrr-\'^\AUcr-')-\ " iMa,L,i < cjr-^''^ A*Jcr-')-' (8) 

and 

\9LrL,i < cgrV-l'^l A:(cr-i)-i, \f\^^^ < cgrV"!'^! exp (-cioA:(cr-i)^) . (9) 

The proof of Theorem 11.31 is an easy consequence of Theorem 11.11 the Hamiltonian ()G2p 
on the domain T>r is, after scaling, equivalent to the Hamiltonian (|G1|) on the domain Di, 
and due to the non-linearity, the small parameter is here naturally r > ^/e instead of e. Note 
that no assumptions on h is required, since no information on the vector 'Vh{L) for / € is 
used except of course at / = 0. 

Let us note that the transformation and the functions g and / are not (a, L)-Gevrey, 
because of the factors r'^^' appearing in the estimates (nevertheless, they are still {a,Lr)- 
Gevrey for some Lr depending on L and r, but this will not be used in the sequel). Even 
though these factors come from our method of proof (which consists in scaling the system so 
that it reduces to a perturbation of a linear integrable Hamiltonian), we believe this is not 
an artefact (similar factors appear in the analytic case when one uses Cauchy estimates to 
control the norm of the derivatives of a function from the norm of the function on a domain 
of size r). 

Let us also note that in the statement of Theorem 11.31 one has the freedom to choose any 
r such that -y/e < r < R, provided r is sufficiently small so that the second part of ([7|) is 
satisfied. In the sequel, we shall use the statement only for r = ^/e which appears to be the 
most natural value, but for some other purposes, for instance if one is interested in global 
stability estimates for non- linear systems (as in the Nekhoroshev theorem), it is useful (and 
perhaps necessary) to use the statement for bigger values of r. 

Now as before, in the Diophantine case we can give a more concrete statement. 

Corollary 1.4. Let LL be as in (jG2p . and uj G 0^(7, r). There exist positive constants c, cq, 
C7, cs, Cg and cio that depend only on n,R,uj,a and L, such that if 

Ve<r <R, r< cCjT^^+^^t, (10) 
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then there exists a symplectic map $ € G"(Pj,/2i ^-r) such that 

Ho<!> = h + g+[f]^ + f, {g,L} = 0. 
Moreover, for L = CL and any I = {li, I2) G N^", we have the estimates 

I'^i - M/L,z,/ ^ cjrr^^^^^ic-'^j-'^r)^, |$e - Ug\^ i i < CTr^l'^l^^-i^-i^^iqi^ 

and 

\s\a,L,i < C8rV~l'2l(c"S"V)iT^, \f\^^^^i < cgrV-'^^l exp (-cio(c7r-i)m) . 

5. It is perhaps worthwhile to compare our results with previous results, which were restricted 
mostly to the analytic case a = 1 or the periodic case d = 1. 

In the analytic case a = 1, results similar to Theorem ll.il (and also to Theorem II .Sp are 
contained in |Pos93j and |DG96j (see also |Sim94j ). Note that our results are however more 
accurate than those contained in [Pos93] : in the estimates for <I> and g, the term A* (c£~^)^^ 
we obtain here is replaced by 1. This does not lead to any improvement if one is interested 
in global stability estimates for perturbations of non-linear integrable Hamiltonians because 
the geometry of resonances prevails in this case, but this improvement on the estimate for 
<I> is visible for global stability estimates for perturbations of linear integrable Hamiltonians 
(compare Theorem 12.11 below with Theorem 5 of [P os93j ) and for local stability estimates for 
perturbations of non-linear integrable Hamiltonians. Also, the improvement on the estimate 
for g is important for the application to the splitting of invariant manifolds (see Remark 13. ip . 
Now concerning Corollary 11.21 in the special case a = 1, one can in fact find better results 
in |Fas90] : still in the estimates for $ and g, the term (c~^7~^e)i+^ we find is replaced by 
c~^7~^e. This discrepancy can be explained as follows. The perturbation theory of quasi- 
periodic motions essentially relies on solving a certain equation (called homological) which 
consists of integrating a function along the linear flow of frequency w. Assuming that the 
function has zero average along the flow, a formal solution always exists, and it is a smooth 
convergent solution if and only if uj is Diophantine. So if one considers non-Diophantine 
vectors (which is important for several problems, such as Nekhoroshev type estimates), one 
can only solve an "approximate" equation. The usual approach is to approximate the function 
by a nicer function (usually a polynomial function), whereas here, as in [BF12j . we just 
approximate the frequency by a nicer frequency (a periodic frequency). But if we restrict to 
Diophantine vectors, then one can actually solve the exact equation, and this ultimately leads 
to better estimates. 

Now in the periodic case d = 1, the result is known for any a > 1 (see [MS0 2]). and our 
proof of the general case 1 < d < n crucially relies on it. In fact, in |MS02j . one has a more 
accurate result, namely the existence of a formal transformation to a formal normal form, both 
having Gevrey asymptotics. Such a result is stronger than the existence of a transformation 
to a normal form with an exponentially small remainder, as the latter is implied by the former 
(but not the contrary). In the analytic case a = 1, these asymptotics were first established 
in [RS96] . 

Finally, there are some related results for d = n and any a > 1. First in [MPlOj . the 
dynamics in the neighbourhood of a Gevrey Lagrangian torus is considered, and a Gevrey 
normal form is constructed. The system studied there can be brought to a system of the 
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form we considered, though with a special perturbation, and the result obtained in [MP10| is 
therefore close but different. Then, in an unpublished manuscript |Sau06j . a result analogous 
to |MS02) . that is the construction of a formal transformation to a formal normal form with 
Gevrey asymptotics, is proved in the case of a Diophantine vector. As before, these two 
results give more information but in a different setting, moreover they are based on solving 
exact homological equations and so they cannot apply to non-Diophantine vectors. 

6. Finally let us describe the plan of this paper. In ^ we will deduce from Theorem 11.11 and 
Theorem [L3] exponentially large stability estimates for the evolution of the action variables, 
which are global for perturbations of linear integrable systems and only local for perturbations 
of non-linear integrable systems. In the first case, that is for perturbations of linear integrable 
systems, we will show on an example that these estimates cannot be improved in general. Then 
in ^ we will use Theorem 11.31 to prove a result of exponential smallness for the splitting of 
invariant manifolds of a hyperbolic tori. The proof of Theorem 11.11 and Theorem 11.31 will be 
given in ^ Finally, we will gather in a appendix some technical estimates concerning Gevrey 
functions that are used to prove our theorems. 

To simplify the notations and improve the readability, when convenient we shall replace 
constants depending on n, R, uj, a and L that can, but need not be, made explicit by a •, that 
is an expression of the form u<v means that there exists a constant c > 0, that depends on 
the above set of parameters, such that u < cv. Similarly, we will use the notations u-> v and 
u = -v. 

2 Application to stability estimates 

In this section, we will give direct consequences of our normal forms Theorem 11.11 and Theo- 
rem [T?3] to the stability of the action variables. 

Recall that we have defined F as the smallest rational subspace containing the vector u. 
The key point to obtain stability estimates is the following observation. By definition of the 
Poisson Bracket, the equalities /cj} = {gJuj} = imply that is a first integral of the 

normal form l^^ + [/]^ + 5 in the statement of Theorem ll.il (or a first integral of /i + [f]i^ + g, 
since h is integrable, in the statement of Theorem ll.3|) . Therefore the action variables of all 
solutions of the normal form stay constant along F, that is Hpilit) — Iq) = as long as 
they are defined, if /q = /(O) is an arbitrary initial action and where we have denoted by 
Hp the projection onto the subspace F. Now if we add a small remainder to the normal 
form, the quantities \IlF{I{t) — Io)\, as long as they are defined, will have small variations 
for an interval of time which is roughly the inverse of the size of (the partial derivative with 
respect to 9 of) the remainder. Coming back to our original system, since it is conjugated by 
a symplectic map which is close to identity to such a normal form with a small remainder, the 
same property remains true as long as the solution stays on the domain where the conjugacy 
is defined. 

1. Let us state precisely the result, starting with a perturbation of a linear integrable 
Hamiltonian as in (|Gip . 

Theorem 2.1. Under the notations and assumptions of Theorem let {9{t),I{t)) be a 
solution of the Hamiltonian system defined by H, with {9q,Iq) G T^r/a o,nd let Tq he the 
smallest t e]0, +oo] such that {9{t),I{t)) ^ 'Dr/a- Then we have the estimates 

\UF{I{t) - Io)\ < 5, \t\< min{ro,fe-iexp(-A:(-e-i)^)} 
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for any (A* (•£ ^))^^ <• S <• 1. Moreover, if F = M", then 

\I{t) - Io\ < S, \t\ <• 6e-^ exp (•A:(-e-i) < 



Corollary 2.2. Under the notations and assumptions of Corollary let {9{t),I{t)) be a 
solution of the Hamiltonian system defined by H, with {0q,Iq) £ T^R/i and let Tq be the 
smallest t £]0,+oo] such that (9{t),I{t)) ^ I'_r/4. Then we have the estimates 

\UF{I{t) - Io)\ <S, |t| <• min|To,(5e"^exp (^•(7e"i)^^^)| 

for any (7" e) 1+^ <• 5 <• 1. Moreover, if F = R", then 

\I{t) - Io\ < S, \t\ <• Se"^ exp (■{-fe~^)^^(^ 



Note that we have exponential stability for the action variables only in the case F = R", 
that is when u is non-resonant, since in this case we can ensure that Tq is at least exponentially 
large. In general, the stability result is only partial since it involves the projection of the 
action components onto F, and the time Tq, which is easily seen to be always at least of order 
6e~^ , cannot be much larger in general (it is easy to construct examples for which the action 
variables of a solution drift along the orthogonal complement of F with a speed of order e). 

Note that in the particular case where 6 = - (A* (-e"^))-!, we have 

\uHm-io)\<{Ki-e'')r' 

for 

|t|<-min{ro,(A:(.e-i))-Viexp(.A:(.£-^)^)} 

1 

and in the Diophantine case, for 6 = ■ (7" e) 1+^ , we have 

\UF{I{t)-Io)\<{j~^e)^ 

for ^ 

\t\ <■ min jro, •7"TT7£-TT7 exp (^■{'J£~^)"^^^^ ■ 

Still in the Diophantine case, we expect that it could be possible to take 6 = -j~^e (see the 
comments at the end of the first section). 

The deduction of Theorem [6] from Theorem 11.11 is very classical, but for completeness we 
give the details. 

Proof. Using Theorem 11.11 it is enough to prove the statement for a solution {6(t),I{t)) of 
the Hamiltonian H = H o ^, with (6*0, Iq) £ ^0 being the smallest t g]0, +00] such 

that {6{t),I{t)) ^ 'Dji/2- Indeed, we have H = H o ^, where $ : 'Dji/2 T^R is symplectic. 
Moreover, from ([5]), we have the estimate 

1^ - Idbo(25^/,) <\^- <■ {AU-e-'))-\ 

so we can ensure that the image of $ contains the domain 2?/j/4, and since we have assumed 
that (A* (-e""^))"^ <• (5 <• 1, the statement for H trivially implies the statement for H (only 
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with different implicit constants). Note that in fact we wiU not use the first part of ([6]), but 
only the following obvious consequence of the second part of ^ : 

\def\cOiv^,,) <■ \f\,L <• eexp (- • iKi-e-'))^) . (11) 

Recall that the equation of motions of H are 

9 = diH{e,i), i = -deH{e,i) 

and using the mean value theorem and the fact that /^j is integrable, 

\I{t)-Io\ < \t\\deH\co(^Vn^,) < \t\\de{[f]u.+g + f)\coiv^^,), |t| < Tq. 

But {g,lu}} = is equivalent to deg{9,I) G F"*-, and similarly for [f]ui, hence if we project 
onto F, we have 

|n^(/(t) - /o)| < \t\\de~f\c^iv^„), \t\ < To. 

Now using ()lip and 

\t\<5e-'exp[iAU-e-'))^) 

we obtain 

\UFim - Io)\ < 6, \t\< min{ro,fe"iexp(-(A:(-e-i))i)} 

for any (A* (•e^^))^-'^ <5< 1. Finally, note that if F = M", the map Hp is the identity so 
that 

ro->fe~iexp(-.(A:(.e~i))i) 

and hence 

\I(.t) - Io\ < 5, \t\ <■ fe-i exp (- • {AU-e-'))^) . 
This concludes the proof. □ 

2. Now we will show on an example that the estimates of Theorem 12. H and therefore the 
estimate on the remainder in Theorem ll.il are essentially optimal. 

Theorem 2.3. Let uj G M"\{0}. Then there exist a sequence {£j)j^fq of positive real numbers 
and a sequence {fj)je_N of functions in G"''^{Vji), with 

£j <■ \fj\a.L <■ £j, lim Ej = 0, 

J— >-+oo 

such that for j •> 1, the Hamiltonian system defined by Hj = luj + fj has solutions {0{t), I{t)) 
which satisfy 

\t\e, exp (- • iAU-ej'))i) <• |n^(/(t) - /o)| <• exp (- • (A:(.eTi))i) . 

Proof. First recall the following fact from linear algebra: there exists an invertible square 
matrix A of size n with integer coefficients such that u = A{'w, 0) G M'^ x = with 
w non-resonant. So replacing H hj H o <^a, where : T" x M" -)> T" x M" is the 

linear symplectic map given by <I>yi(0,/) = {AO, (A~^)*/), we can assume that w = (ti7,0) as 
the statement for H and H o cj)^ are equivalent (up to constants depending on A, and so on 
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n and u). Moreover, dividing w by ±|ci7| and reordering its components if necessary, we can 
also assume that w = (1, ro, . . . , Wd-i), with \zui\ < 1 for 1 < i < d — 1. 
So we have reduced the general case to 

oj = (l,zi7i,... ,CTd-i,0) G M'' X M"" = M", F = M.'^ X {0}. 

Let us denote by {pj/qj)j£n the sequence of the convergents of zui for instance. We have the 
classical inequalities 

(Qj + Qj+i)''^ < ki^i -Pj\ < Qj+v 3 e N, 
and since ^j+i > qj, this gives 

{2qj+i)-^ <\qjWi-pj\<q-l^, j E N. (12) 
Now by definition of we obtain 

qj+l < ^^{qj) < 2qj+i. (13) 
The perturbation fj will be of the form 

fjie,i) = f}{i) + ff{e), ie,i)GVR. 

First, we choose fj{I) = Vj ■ I — uj ■ I, where vj = {l,pj /qj, zu2, . . . , Wd-i,0). We set 

Ej = \wi -Pj/qj\. 

From the inequalities (fT2]) and (fT3]l . recalling the definitions of A^^ and A*, we have 

(A^(<?,-))-' <e, <• (A^(g,))-\ a:(s^') < '^.^ < a:(-s"')- (14) 

Now by definition, \fj\a,L=- £j- Then, if we let kj = {pj, —qj,0, . . . ,0) G Z", we define 
fj{0) = EjHj sm{27rkj ■ 6) with ^ij to be chosen. Note that \kj\ = qj since \qj\ > \pj\ (as 
1^1 1 ^1)- It is easy to estimate 

Sjfj.jexpi-lkjl^/") <■ \fj\a,L <-ejfj.jexp{-\kj\^/°') 

and so 

ejfijexp{-qj' ) < \fj\a,L < exp{-qj' ). 

Therefore we choose fij =■ exp(— • Q'j''") so that £j <• l/J |a,L <"^j- Finally, fj = fj + /? so 
£j <• <• Ej, and — )• when j — )• +oo. Now we can write the Hamiltonian 

Hj{e,i) = L{i) + fj{e,i) 

= uj-I-\-Vj-I — uj-I-\- £jfij sin(27r/!:j • 9) 
= Vj ■ I + Sjfij sin(27rA;j • 9) 

and as kj.Vj = 0, the associated system is easily integrated: 

r 9{t) = 9o + tVj [Z"] 

\ I{t) = Iq — t2-7Tkjejfj,j cos{2Trkj.9(j). 
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Choosing any solution with initial condition (9q,Iq) satisfying kj.Oo = 0, cos{2Trkj .Oq) = 1 
and using the fact that \kj\ = qj and /Xj =• exp(— • q^"'), we obtain 

IHt) - ^ol =• \t\£jqj exp(- • qy^). 

For j •> 1, Qj •> 1 and we have exp(— • gj''") < qj exp(— • (/j^") < exp(— • ^j^") for well-chosen 
implicit constants, hence 

\t\ej exp(- • qy^) <• \I{t) - Io\ <• \t\ej exp(- • gj^"). 

Using ([HD this gives 

\t\e, exp (- • {A*J.e-'))^) <■ \I{t) - Io\ <• \t\e, exp (- • {A*J.e-'))^) 

and as the vector kj belongs to F = M'^ x {0}, the solution we have constructed satisfy 
UF{I(t) - Jo) = I{t) - Iq. This concludes the proof. □ 

3. For a perturbation of a non-linear integrable Hamiltonian as in ()G2p . we have a stability 
result similar to Theorem 12.11 

Theorem 2.4. Under the notations and assumptions of Theorem \1.3\ let {6(t),I{t)) be a 
solution of the Hamiltonian system defined by H, with {9q,Iq) £ T^r/A ^^'^ ^ei Tq be the 
smallest t g]0,-|-oo] such that {9{t),I{t)) ^ 'Dr/4- Then we have the estimates 

\UFim -Io)\< S, \t\ <• min {To, 6r~^ exp (•A:(t-1)^) } 

for any r(A* (-r^^))^^ <• 6 <• r. Moreover, if F = R", then 

|/(t) -Io\<6, \t\ <• Sr-^ exp (^-AU-r-^)^) ■ 

Corollary 2.5. Under the notations and assumptions of Corollary \1.4[ let {9{t),I{t)) be 
a solution of the Hamiltonian system defined by H, with {6q,Iq) £ 'Dr/A ^^'^ be the 

smallest t g]0, -l-oo] such that {9{t),I{t)) ^ Then we have the estimates 

|nF(/(t) - /o)| <• S, \t\ <• min |ro, 6r-^ exp (^■{-fr-^)^^(^^^ 

for any r(7~ r) <• 6 <• r. Moreover, if F = M", then 

\I{t) - Io\ <• 6, \t\ <• 5r~2 exp (^■{jr^^)^^(^'^ . 

The proof of Theorem [23] is entirely similar to the proof of Theorem [57T1 replacing the use 
of Theorem II. II bv the use of Theorem 1 1.31 since the only information we need to derive these 
estimates from the resonant normal form with a remainder are estimates for the distance 
to the identity (with respect to the action variables) of the conjugacy $ and estimates 
for the partial derivatives (with respect to the angle variables) of the remainder /, and these 
information are contained in Theorem 11.31 Therefore we shall not repeat the details. 

These local estimates of stability (or variants of them) are at the basis of global (in phase 
space) estimates of stability for perturbations of non-linear integrable systems, provided that 
the integrable part satisfy some geometric assumption. However, we don't know if it is 
possible to construct an example such as in Theorem 12.31 that would show that the estimates 
of Theorem 12.41 are "essentially" optimal. 
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3 Application to the splitting of invariant manifolds 



In this section, we apply our normal form results to a different but ultimately related problem, 
which is the so-called "splitting" of invariant manifolds (or "splitting" of "separatrices" ) . 
Roughly speaking, if a Hamiltonian system H as in ()G2p has a suitable "hyperbolic" invariant 
torus for which the stable and unstable manifolds intersect, the problem is to evaluate in 
some sense the "angle" between these invariant manifolds at the intersection point. This is 
an important problem in itself, and this is also deeply connected to the problem of the speed 
of instability for the action variables (which is known as the speed of Arnold diffusion). 

The general principle is that the "splitting" is exponentially small for analytic systems 
and the literature on the subject is huge. Here we shall closely follows [ LMS03J , Chapter §2, 
where an approach to obtain exponentially small upper bounds in the analytic case is given 
based on normal forms techniques. The results contained in |LMS03] are quite general, as 
they are valid for any number of degrees of freedom, any Diophantine frequency vector and 
without any restriction on the perturbation, assuming the torus exists and that its stable and 
unstable manifolds intersect. However, because of this great generality, these results are not 
very accurate, they only give an upper bound with a reasonable value for the exponent in the 
exponential factor. Much more accurate results, such as an asymptotic formula and reasonable 
values for the other constants involved, can be obtained in much more restricted situations 
(essentially for two degrees of freedom, specific frequency vectors and specific perturbations, 
see [BFGS12] for some recent results and references). 

Our aim here is to generalize the results of [LMS03j for Hamiltonian systems which are 
only Gevrey regular. We will also assume the existence of a "hyperbolic" torus, together 
with the property that its invariant manifolds intersect (that is, the existence of homoclinic 
orbits), simply because conditions that ensure the existence of these objects are well-known 
(we will quote some of these conditions). 

1. We will not try to give an abstract definition of a hyperbolic torus for a Hamiltonian 
system, first because we will deal with rather concrete examples, but also because we could 
not find any satisfactory abstract definition (for instance one which would ensure the existence 
and uniqueness of stable and unstable manifolds, we refer to |LMS03| and |BTOO) for some 
attempts). 

Now let us consider the setting as described in ()G2p . Let To = x {0} be the invariant 
torus, for the integrable system, with frequency oj G M" \ {0}. Without loss of generality 
(as we explained in the proof of Theorem 12. 3|) . we may already assume that oo = (ti7,0) G 
R-i X M™ = R" with w non-resonant. 

If a; S r2„(7,r), by KAM theory this torus will persist under any sufficiently small and 
regular perturbation, provided that h is non-degenerate (in the analytic case, the assumption 
that OJ G l^n(7j ''") can be weakened). This is not true \l d <n — l; however, under appropriate 
assumptions on the system, one can still apply KAM theory to prove the existence of an 
invariant d-dimensional torus, with frequency zu, which is hyperbolic in the sense that it 
possesses stable and unstable manifolds. Moreover, such a tori will be isotropic and its 
asymptotic manifolds will be Lagrangian. If the stable and unstable manifolds intersect, 
following [LMS03] . we can define a symmetric matrix of size n at a given homoclinic point, 
called a splitting matrix, the eigenvalues of which are called the splitting angles. Our result is 
that there exists at least d splitting angles which are exponentially small (see Theorem 13. 3|) . 
The proof of this result will be analogous to the one in |LMS03j . However, in [LMS03] . it 
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seems that the normal forms they used, which are taken from |Pos93] . are not accurate enough 
to derive the result they claimed (see Remark 13. ip . 

Before stating precisely the result, we need some rather lengthy preparations in order 
to introduce our definitions and assumptions (the validity of our assumptions will be briefly 
discussed after the statement of Theorem I3.3p . 

First let us split our angle- action coordinates (0,1) S Vr = x Bji accordingly to 
w = (ro, 0) e M*^ X = M": we write d = (^i, ^2) G x and / = (Ii, /z) e Bf^ x 5^, 
where 3^^ = BrH and = 5r n M*". We shall always assume here that m > 1, that is 
d < n — 1. Then Theorem 11.31 states that for the value r = 2i/e, if 

there exists a symplectic map $ € G"(P^/2 5^r); satisfying the estimate dSD, such that 
Ho<^ = h+[f]^ + g + f, [fUe2,I)= [ f{9,I)dei, g{9J)=g{62j) 

with g and / satisfying the estimates dH). As in the proof of Theorem \1.3\ it will be more 
convenient to use a rescaled version of the Hamiltonian H o ^. Consider the map 

a: {9,1) ^{9,^1) 

which sends the domain Di onto T>j./2, and let 

n = £'^{H o^oa) = e~^h o a + £'^[f]u o cr + e'^g o cr + e'^ f o a 

be the rescaled Hamiltonian, which is defined on Note that in the proof of Theorem 11.31 
we used the same scaling map but the scaling factor was y/e ^ in order to bring the linear 
part of h to order one, while here the scaling factor is since it is the quadratic part of h we 
want to bring to order one. Consequently, the solutions of the Hamiltonian system defined by 
% coincide with those of H only after scaling time by y/e: since we will be only interested 
in invariant and asymptotic manifolds, it makes no difference to consider % instead oi H o^. 
Now let us define 

A = A(e) = (A:(.V^-'))~\ /i = ^(e) = exp(- • A(e)-i) = exp(- • A"^). (15) 
For any I = (hjh) ^ N^", we have 

d^{g o a) = d'id\\g o a) = V^^''\d'g) o a 

so 

|5'(5°f^)lco(Bi) < \/^''''|5'5'lco(C,/2) 

and therefore 

Now using Q, this gives 

l5°^L,L,z<•^^"''e(^/^)""''A = •eA 
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and hence 

\9°(^\a,L<-^^- (16) 

In the same way, 

\foa\^j^<eti. (17) 

We win further decompose h and [fj^j as follows. Without loss of generality, we can assume 
that /i(0) = 0. Then by Taylor's formula, we can expand /i at / = at order 2: 

h{^/eI) = y/ezu ■ h + eAh ■ h + eBh ■ h + eCh ■ h + ey/eRh{I) 

where 



Rh{I) = 2-^ [\l-tfV^h{tV^I){I,I,I)dt, 
Jo 



A and B are square matrix of size respectively d and m, and C is a matrix of size m times d. 
The rectangular term represented by CIi ■ I2 is slightly inconvenient, so we will assume: 

(A.l) C = 0. 

Then we expand [f]^^ at / = at order 0: 

[fU02, V^i) = [fU02,o) + ^ C V/[/]^(e2, t^i)idt = ¥{62) + (02, /). 

^0 



and so 

Thus we have 



n{e,i) = y/l w-h + Ah-h + Bh-h + v{e2) 

+ {V~eRh{I) + J02, /) + £-'5(^2, V~el)) + e-'fie, ^el). 

Applying Lemma IA.3I (respectively with s = 3 and s = 1), we obtain \Rh\a,L' <■ 1 and 
\R[f]^\a.L' <^ ^1 for L' = L/2, so that the first two terms in the parenthesis above (which 
are the Hamiltonians independent of ^i) are of order ^/£. However, the last term in this 
parenthesis is of order A > ^/£ by (jl6p . so we define 

R = \-^^R^ + \-^^'^R^f^^ + \-^£-^g o CT, 

so that the Hamiltonian in the parenthesis above is Ai?, with <• 1. Let us also write 

F = e'^^i~^foa 

so that the last term in the expression of % is ^F, with < <• 1 by (fT7|) . Finally, 

we obviously have |T^|a,L' < 1- 

Now recalling the dependence of A and ^ in the notation, we have obtained a Hamiltonian 
% = ^A,/i5 defined on of the form 

f UxA'^, I) = nx{e2,i) + fiF{e, /), \f\^^l' <• i 
nx{e2j) = w{e2j) + \R{e2j). \R\a,u<i (is) 
[7^-(02, /) = V^'^^ ■ h + Ah ■ h + BI2 ■ I2 + v{92), \v\^,L' < 1. 
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It will be convenient to consider ^a,^ as a Hamiltonian depending on two independent pa- 
rameters A and fi, and as a rule the notation A(e) and ^{e) will be used only when we want 
to recall that they both depend on e as in (jlSp . One has to consider Tix^/^ as an arbitrary 
//-perturbation of Tlx = Hxfl, and Hx as a special A-perturbation of the "averaged" system 
Ti""'" (special because the perturbation R is independent of ^i). Note that the averaged system 
can be further decomposed as a sum of two Hamiltonians 

n''''{92,I) = K{h) + P{92,l2), 

where K{Ii) = ^/e ■ I\ + AI\ • /i is a completely integrable system on = T"' x B'^ and 
P{02,l2) = BI2 ■ I2 + ^(^2) is a mechanical system (or a "multidimensional pendulum") on 
= T™ X We make the following assumption on the mechanical system: 

(A. 2) The matrix B is positive definite (or negative definite), and the function V : T'" — t- M 
has a non-degenerated maximum (or minimum). 

Without loss of generality, we may assume that V reaches its maximum at ^2 = 0, so that 
O = (0, 0) G Pj" is a hyperbolic fixed point for the Hamiltonian flow generated by P. This in 
turns implies that, for any /* E Bf, the set T{I*) = {/i = Ii} x O is a d-dimensional torus 
invariant for the averaged system, and it is hyperbolic in the sense that it has stable and 
unstable manifolds 

w^inn) = {h = 1*1} X w^{o) 

where W^{0) are the stable and unstable manifolds of O, which are Lagrangian. In particular, 
the torus T(0) is quasi-periodic with frequency y/e ^w. 

Now this picture is easily seen to persist if we move from T-L""'" to T-Lx- Indeed, since T-Lx is 
still independent of 0i, the level sets of Ii are still invariant, hence for a given G Bf, the 
Hamiltonian flow generated by the restriction of Hx to {Ii = If} x I?™ (considered as a flow 
on depending on /^) is a A-perturbation of the Hamiltonian flow generated by P: as a 
consequence it has a hyperbolic fixed point Ox{Ii) G which is A-close to O, for A small 
enough. Hence Tx{I*) = {Ii = Ii} x Ox{Ii) is invariant under the Hamiltonian flow of Hx, 
and it is hyperbolic with Lagrangian stable and unstable manifolds 

w^iTxiH) = {h = m X W^iOxilD). 

The torus 7a(0) is still quasi-periodic with frequency ^/£ ^w. 

Remark 3.1. Let us point out here that in ILMSOSij . it is stated incorrectly that the size of 
the Hamiltonian we called \R is of order (as we mentioned above, this Hamiltonian is 
composed of three terms and only two of which are of order y/e, the last one being of order 
X{e)). Now in JLMSOS^ . they made use of normal forms taken in JPds93^ which contains 
estimates that are less accurate than ours (as we already explained, A(e) is only of order 
one in JPds93^ ) and consequently these estimates do not allow to show the existence of the 
invariant torus for the system we called Hx ■ 

Our next assumption concerns the persistence of the torus 7a (0) , as well as its stable and 
unstable manifolds, when we move from Hx to Hx,f_i'- 

(A. 3) For any < A <• 1 and < ^ <• A, the system ^a,^ has an invariant torus 7a,^, with 
7a,o = Tx = TxiO), of frequency ^/e ^w, with stable and unstable manifolds W^{Tx,fx) 
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which are exact Lagrangian graphs over fixed relatively compact domains U C T". Moreover, 
W'^{Tx^ij,) are /x-close to W'^{T\) for the C^-topology. 

Let us denote by Sf^ generating functions for W^(Tx^^) over U^, that is if = x Bi, 
then 

where : — ?■ M are functions. Since W^{Tx^^) are /i-close to W^{Tx) for the 
C^-topology, the first derivatives of the functions S^^ are //-close to the first derivatives of 
Sf = S^Q for the C^-topology. 

Then, in order to evaluate the splitting, we need the existence of orbits which are homo- 
clinic to 7a,^: 

(A.4) For any < A <• 1 and < ^ <• A, the set VF+(7^,^) n W-{Tx,^,) \ Tx,^ is non-empty. 

2. We can finally define the notions of splitting matrix and splitting angles, and state our 
results. 

The set {Tx,n)f~\W^ {J'x,ti)\Tx,^ is invariant, so it consists of orbits of the Hamiltonian 
system defined by 'Ha,/^. Let 7^,^ be one of this homoclinic orbit, and px,ii = 7a,^(0) = 
{9x,fi, Ix,iJ,) a homoclinic point. Since px,^ is a homoclinic point, 9x,n G n and 
^e5'^^(6'A,Ai) = deS^^^{9x,^)- Then we can define the splitting matrix M{Tx,^,Px,^J.) of Tx,^ 
at the point px,fj,, as the symmetric square matrix of size n 

M(rA,^,PA,M) = dlisl^ - 5,-^)(^a,m). 

Moreover, for 1 < i < n, we define the splitting angles ai{Tx,fj.,Px,n) as the eigenvalues of the 
matrix M{Tx,fi,Px,fi)- Note that since the homoclinic point belongs to a homoclinic orbit, at 
least one of these angles is necessarily zero. 

Theorem 3.2. Let 'Hx,^i be as in p^ . and assume that (A. 2), (A. 3) and (A.4) are satisfied. 
Then, with the previous notations, we have the estimates 

\ai{Tx,^i,P\,tJi)\< ^J^, l<i<d. 

We have used the fact that the stable and unstable manifolds are exact Lagrangian graphs 
over some domains in T" in order to define the splitting matrix and splitting angles, but in 
fact only the Lagrangian property (and not the exactness nor the graph property) is necessary 
to make those definitions (see [LMS03J). 

Now the solutions of the Hamiltonian system defined by if o $ differs from those of 
'Hx(e),i_i{e) only by a time change, so Tx(e),^{e) is still an invariant hyperbolic torus for H o 
with the same stable and unstable manifolds. Coming back to our original system, the 
torus Te = ^{Tx{e), 11(e)) is hyperbolic for H, with stable and unstable manifolds W^{Ti;) = 
^{W^iTx{e),ii{e))), and for 7^ = $(7A(e),Ai(e)) and Pe = $(PA(e),/.(£)) we can define a splitting 
matrix M{Te,Pe) and splitting angles ai{Te,Pe) for 1 < i < n. 

Theorem 3.3. Let H be as in (|G2p . with r = satisfying ([7]). Assume that (A.l) is 
satisfied, and that (A. 2), (A. 3) and (A.4) are satisfied for the Hamiltonian 'Hx{e),ii{e)- Then, 
with the previous notations, we have the estimates 

\ai{Te^Pe)\ <• (1 + (A:(.V^-'))-i) exp (- • A:(.V^-')^) , l<i<d. 
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Corollary 3.4. Let H he as in ()G2p . withuj G 1^(^(7, r) andr = 2y/£ satisfying (jlOp . Assume 
that (A.l) is satisfied, and that (A. 2), (A. 3) and (A. 4) are satisfied for the Hamiltonian 
^A(e),/i(e)- Then, with the previous notations, we have the estimates 

\ai{Te,Ps)\ <V^{l + (-T-^e)^^) exp (- • (^e"^)^^^) , I < i < d. 

Let us note that the exponent in the exponential factor in our corollary is not far from 
being optimal, at least for a>l, d = n — 1 and for a badly approximable vector oj £ 
r2„_i(7, n — 2): indeed, in this case we have an exponentially small upper bound with the 
exponent (2a(n — whereas in |Mar05) . a sequence of ej-perturbations, with ej going 

to zero as j goes to infinity, is constructed such that the perturbed system has an invariant 
hyperbolic torus of dimension d = n — 1, with d — 1 = n — 2 splitting angles which have an 
exponential small lower bound with the exponent (2(a — l)(n — 2))^^. 

Let us now briefly discuss the validity of our assumptions (A.l), (A. 2), (A. 3) and (A. 4), 
referring to |LMS03j for more details. Concerning (A.l), in principle it is just a simplifying 
assumption and it can be removed, though we shall not try to do it here. The assumption 
(A. 2) is crucial as it ensures that the averaged system has an invariant hyperbolic tori. Then, 
using classical KAM theory, (A. 3) follows from (A. 2) under usual assumptions (there are 
many references in the analytic case, but we do not know any for non-analytic but sufficiently 
regular systems). Now concerning the existence of homoclinic orbits (A. 4), this is a general 
assumption as follows: using a variational argument one can prove that the mechanical system 
has orbits homoclinic to the hyperbolic fixed point and so the averaged system has orbits 
homoclinic to the hyperbolic torus 7^(0), then assuming that the stable and unstable manifolds 
of T(0) intersect transversely along one of this homoclinic orbit inside the energy level, the 
assumption (A. 4) is satisfied, that is this homoclinic orbit for the averaged system can be 
continued to a homoclinic orbit for the full system. 

Note that the results in Theorem l3.2l and Theorem [331 are just upper bounds on d splitting 
angles. These splitting angles can be actually equal to zero, as we did not assume that the 
stable and unstable manifolds intersect transversely inside the energy level. As we already 
mentioned, this transversality assumption implies (A. 4) provided (A.l), (A. 2) and (A. 3) are 
satisfied, and moreover, under this assumption, one can prove that they are exactly d splitting 
angles which are non-zero and exponentially small as the other angles can only be polynomially 
small (we recall that at least one of them is zero, but we could have avoided this situation by 
taking a Poincare section and studied the associated discrete system). 

3. Let us now give some details concerning the proof of Theorem 13. 21 and Theorem l3.31 which 
follows from simple lemmas (the proof of which, or references, can be found in |LMS03| ). 

Consider first the Hamiltonian Hx, with its hyperbolic invariant torus 7a = 7a, = 7a (0) 
with stable and unstable manifolds 

W^iTx) = {h = 0} X W^iOxiO)). 

Let us denote by Sf = Sf^ the corresponding generating functions, and by px = Px,o the 
homoclinic point. From the expression of W^{Tx), one immediately has de-^^S^ = 0. Let us 
write Ma = M(Tx,Px)- The following lemma is then obvious. 

Lemma 3.5. The matrix Mx admits the following block decomposition: 

/ \ 



18 



where = dg^{S'^ — )(^a) is a square symmetric matrix of size m. 

Note that the matrix always have zero as an eigenvalue (because we are in a continuous 
setting), and that the fact that W^{Tx) intersect transversely along the homoclinic orbit inside 
the energy level is equivalent to zero being a simple eigenvalue of M^. 

Now let us come back to the Hamiltonian H\ fj^, and let us write Mx^^ — Af(7A,^,PA,/i)- 
From assumption (A. 3) the second derivatives of the functions S^^ are /i-close to the second 
derivatives of for the C'^-topology. This immediately implies the following lemma, where 
we also denote by | . | the norm induced by the supremum norm on the space of matrices. 

Lemma 3.6. We have the estimate 

\Mx,i, - Mx\ < fi. 

Then, in order to use the previous lemma, we need to know how the eigenvalues of a 
symmetric matrix vary under perturbation. This is the content of the next lemma, where we 
denote by d the distance on induced by the supremum norm. 

Lemma 3.7. Let A and A' be two symmetric matrices, with spectrum Spec(^) and Spec(^'). 
Then we have the estimate 

(i(Spec(A),Spec(^')) < l^-^'l- 

Let us remark that the above lemma is quite specific to symmetric matrices, and the fact 
that the splitting matrix is indeed symmetric in the general case (in our restricted case, this 
is obvious by definition) ultimately comes from the Lagrangian character of the stable and 
unstable manifolds. 

The proof of Theorem [32] is now a straightforward consequence of Lemma [331 Lemma [321 
and Lemma 13.61 

Let us now come back to the Hamiltonian H as in (jG2p . It is related to the Hamiltonian 
^A(£),/i(e) by a scaling transformation and by a symplectic transformation <I>. The effect of 
the scaling transformation is simply to multiply the splitting matrix by ^/£. The effect of the 
symplectic transformation is more complicated to describe, but the overall effect is contained 
in the following lemma, where we denote = M{Ts,Pe)- 

Lemma 3.8. There exist two invertible square matrix B and C of size n such that 

Me = V^BMx^e)Me)C 

where B and C satisfy 

\B\ <• 1 + (a:(.^-^))-i, \c-^\ <• 1 + (a:(.v^-^))-i. 

The matrices B and C depends on the differential of the transformation $ (we refer to 
|LMS03j . Proposition 1.5.4 for an explicit expression), and the estimates on B and C follows 
from the estimates on de^Q and dg^i contained in ([8]) in Theorem 11.31 

Finally, we need yet another lemma from linear algebra. 

Lemma 3.9. Let A and A' two symmetric matrices, with A' = BAC for two invertible square 
matrix B and C of size n. Assume that A has d eigenvalues in the interval {—5,5). Then A' 
has d eigenvalues in the interval {—5\B\\C~^\, 5\B\\C~^\) . 

The proof of Theorem 13.31 is now an obvious consequence of Theorem 13.21 Lemma [3. 8 1 and 
Lemma 13.91 
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4 Proof of the main results 



This section is devoted to the proofs of Theorem 11.11 and Theorem II. 3[ 

1. Let us first state an approximation lemma of an arbitrary vector by hnearly independent 
periodic vectors, which was proved in |BF12j . Recah that a vector v E M" \ {0} is cahed 
periodic if there exists a real number t > such that tv G Z". In this case, T = inf{t > 
\ tv & Z"} is called the period of v, and a periodic vector with period T will be simply 
called T-periodic. It is easy to see that a vector is periodic if and only if its minimal rational 
subspace is one-dimensional. 

Now consider an arbitrary vector uj G \ {0}, and let d be the dimension of its minimal 
rational subspace F = F^. For a given Q > 1, it is always possible to find a T-periodic vector 
V £ F, which is a Q-approximation in the sense that \Tu! — Tv\ < Q~^, and for which the 
period T satisfies the upper bound T<-Q'^~^: this is essentially the content of Dirichlet's 
theorem. Then it is not hard to see that there exist not only one, but d linearly independent 
periodic vectors in F which are Q-approximations. Moreover, one can obtain not only linearly 
independent vectors, but periodic vectors wi, . . . , of periods Ti, . . . ,Td such that the integer 
vectors Tivi, . . . , T^Vd form a Z-basis of Z" n F. However, the upper bound on the associated 
periods Ti, . . . is necessarily bigger than Q'^~^, and is given by a function that we call 
here (once again, see [BF12] for more precise and general information, but note that there 

was denoted by ^'^j and ^^j, which we defined in ([1]), was denoted by ^'^). The main 
Diophantine result of |BF12] is that this function is in fact equivalent to the function ^^j, 
up to constants and for Q large enough. This gives the following result. 

Proposition 4.1. Let lo G M" \ {0}. For any Q •> 1, there exist d periodic vectors vi, . . . , v^, 
of periods Ti, . . . , Td, such that Tivi, . . . , TdVd form a Z-basis o/Z"nF and for j £ {1, . . . ,d}, 



\UJ 



v,\<iTjQ)-\ 1<-T,<-^^(Q). 



For the proof, we refer to [BF12j . Proposition 2.3. The implicit constants depends only 
on d and u (the dependence on uj is through its norm \uj\ and the discriminant of the lattice 

Z"nF). 

Now a consequence of the fact that the vectors Tivi, . . . , TdVd form a Z-basis of Z" n F is 
contained in the following corollary. For simplicity, we shall write [ • = ['■■[■ ]fi ]i)d' 

where [-J^ has been defined for an arbitrary vector u) in ([3]). 

Corollary 4.2. Under the assumptions of Proposition ]^. 1\ let = oj ■ I and ly. (I) = Vj ■ I 
for j e {1, . . . ,d}. For any g £ C^{Vr), we have [g]^ = [g]vi,...,Vi and therefore {g^L} = if 
and only if {g, l^. } = for any j £ {1,... ,d}. 

Proof. Let A = Z" n F and dt) be the Haar measure on the compact quotient group F/A. 
The flow leaves invariant the foliation on T" induced by the trivial affine foliation on M" 
defined by F. Each leaf of this foliation on T" is diffeomorphic to F/A, and the restriction of 
to each leaf is uniquely ergodic. By Birkhoff's ergodic theorem, we have 



[g]^ = hm - / 5 o Xidt = [ go X^d^. 
Also, as TjVj is an integer vector, we have 



dt. 
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Using the fact that the vectors Tivi, . . . , T^Vd form a Z-basis of A, the first assertion follows 
easily from these expressions by a change of variables. Now if {g,lvj} = for any j G 
{1, . . . , d}, then g = [g]^^ and so by the first assertion, g = [g\u)-, which means that {g, 1^)} = 0. 
Conversely, if {g.,lui] = 0, then g = [5]^ and therefore g = [g\vi,...,va by the first assertion. 
Since the maps [■]vj are projectors (that is [-Ivj.ijj = ['Ivj), this implies that 

\9\vd — [g]vi,...,vi,vi = \g]vi,...,vi = g, 

and since they commute, we eventually find [gj^^ = g for any j G {1, . . . ,d} and therefore 
{g, Ivj } = for any j G {1, . . . , 4. □ 

2. Then we shall make use of the statement of Theorem ll.il (or Corollarv ll.2p in the particular 
case d = 1, that is when F is one-dimensional, which was proved in [MS02j . As we already 
said, in this situation the vector is in fact periodic so we shall denote it by f , and for any 
non-zero integer vector k £ F, we have the lower bound |A: • u| > and so we will have 
A*(e~-'^) > (Te)"^ in the statement of Theorem 11.11 (or r = and 7 = T"^ in the statement 
of Corollarv ll.2p for this particular case. 

For subsequent use, we introduce another parameter > and we consider the Hamilto- 
nian 

H{9J) = k{I) + s{I) + u{e,I), (9,1) eVR, ^^g^ 
G Z", |s|a,L < z^, \u\a,L<£:- 

Let us define 6 = {2d)^^R. Then we have the following result. 

Proposition 4.3. Let H be as in (|19p . and assume that 

e<iy, Tu<l. (20) 

Then there exists a symplectic map <I> G C^'^' {Vji_s,T^r), with C = -l and V = CL, such 
that 

H0^ = l^ + S+ [U]y +U' +U, {U', ly} = 

with the estimates 

|$-IdU,L/<-rz^, \u'\o,,L' <£Tu, \u\a,L' <£exp{- ■ {Tuy'^/'^). 

For the Hamiltonian (fT9]l . the term is considered as unperturbed, and s -|- w is the 
perturbation. Since we have assumed that e <■ fi, the size of the perturbation is of order n, 
but as s is integrable, the size of the non-resonant part of the perturbation is of order e. Note 
that if we are only interested in the periodic case, then one may take s = in ()19p . 8 = 1^ 
and write u = f in the statement of Proposition 14.31 and this gives exactly the statement of 
Theorem 11.11 

For a proof of Proposition \A.3\ we refer to |MS02j . Proposition 3.2 (note that in |MS02] . 
the size of the perturbation is called e instead of /i, and the size of the non-resonant part 
of the perturbation is called e' instead of e). The implicit constants in the above statement 
depends only on n, it!, a and L. 

Let us now state a simple algebraic property (in our restricted setting), which complements 
Proposition 14.31 
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Proposition 4.4. Under the assumptions of Proposition \4-3\ suppose that 1^ is a linear 
integrable Hamiltonian such that {u,/^} = 0. Then, in the conclusions of Proposition we 
have {u' , lyj] = Q. 

This property has been used several times (this was first used in |Bam99) and |Pos99) ) and 
is valid under much more general assumptions, we refer to [BN12] for a proof in the analytic 
case (see also jBoullj ). Of course, this is a purely algebraic property, and is valid regardless 
of the regularity of the system. 

3. Now we can finally prove Theorem ll.H which is a straightforward consequence of the 
more flexible proposition below. 

Proposition 4.5. Let H be as in ()Gip . and Q > 1. If 

Q->1, e<AUQr\ (21) 
then there exists a symplectic map $ € G°'^^{Vjii2,Vr), with L = -L, such that 

Ho<^ = l^ + [f]^ + g + f, {g,L} = 

with the estimates 

\<l>-ldlj^<Q-\ bL,z<-eQ-\ l/UL<-^exp(-.Q-V«). 
Proof of Theorem We choose 

with a well-chosen implicit constant so that the second part of (|2ip is satisfied. Proposi- 
tion 14.51 with this value of Q implies Theorem II. H as the first part of ()2ip is satisfied by the 
threshold (jH). □ 



It remains to prove Proposition 14.51 

Proof of Proposition Recall that we are considering H as in (|Gip . Since Q ■> 1 by the 
first part of ([2T]) . we can apply Proposition 14.11 there exist d periodic vectors vi, . . . ,Vd, of 
periods Ti, . . . , T^, such that Tivi, . . . , TdVd form a Z-basis of n F and for j G {1, . . . , d}, 

\uj - Vj\ <■ iTjQy\ 1 <• Tj <• ^UQ)- 

For j £ {1, ... ,d}, let us define 

Sj — luj ^Vj ) l^j — ■ (I^jQ) 

with a suitable implicit constant so that |sj|a,L < t^j- Note that l^j = Iv- + sj, and that 
Tji^j =■ Q~^. Let us further define Lj = C^L and Rj = R — j6 so that in particular R^ = R/2. 

Then we claim that for all 1 < j < d, there exists a symplectic map S C^'^^ {T>r.,T>ji) 
such that 

Ho<^j = l^ + [f]y,^,„^y^ + gj + fj, {gj,lvj = 
for any 1 < i < j, with the estimates 

\^j - Id|a,L^. <Q~\ \gj\a,L, < sQ'^ , <-eexp(- • Q-^/°). 
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The proof of the proposition fohows from this claim: it is sufficient to let ^ = ^>^, g = g^, 
f = fd and L = L^, since from Corollary [f]vi,...,Vd = [f]uj and {g, /„J = for 1 < z < d 
is equivalent to {g, = 0. 

Now let us to prove the claim by induction. For j = 1, writing l^^ = /^^ +si, this is nothing 
but Proposition 14.31 (up to a change of notations, namely si instead of s, f instead of u, gi 
instead of u' and /i instead of u). So assume the statement holds true for some 1 < j < d — 1, 
and let us prove it is true for j + 1. By the inductive assumption, there exists a symplectic 
map G C^'^i [Vr^ , Vr) such that 

H o^j = l^ + [f]v^,...,v, + gj + fj, {gj,lv,} = 0, 

for any 1 < » < j, with the estimates 

l^i - U\a,L, <Q^\ \9j\a,L, < eQ"\ \fj\a,L, < E exp{- ■ Q~^/°). 

Now let us consider the Hamiltonian 

H o - fj = + [/]t,i, + 9j = Ivj+i + Sj+l + [f]v^,...,Vj + gj = kj+j + Sj+i + Uj 

with Uj = [f]vi,...,vj +gji so \uj\a,Lj < {£+sQ~^) < £■ We want to apply Proposition l4.3l to this 
Hamiltonian, and we observe that ([21]) implies ([20]) : indeed, Q •> 1 implies TjVj =■ <• 1, 
whereas 

e <• AUQ)-' =■ (Q'fUQ))'' <• {TjQ)-'=-Uj. 

So we can apply Proposition 14.31 there exists a symplectic map <I)-?+^ £ C^.^j+i (T>r.^^ , Vr. ) 
such that 

{H o $j - fj) o = + sj+i + [uj]^^^, + u'j + Uj 

= L+ [[f]vi,...,vj + 9j]vj+i + u'j + Uj 

= L + [f]vu-,vj+i + [gj]v,+i + u'j + Uj 

with {u'j, Ivj+i} = and the estimates 
and 

\u'j\a,Lj+i <£Q~^, \uj\a,Lj+i <-eexp(- • (Tj+ii/j+i)"^/") =-eexp(- • Q"^/"). 
Obviously, we have | [(7j]t;j_,_i |a,Lj <£Q~^, so we set 

$^.^1 = o gj+i = [gj]^^^^+u'j, fj+i=Uj + fjO^^+\ 

We have ^j+i £ 0°''^^+^ {Vr^^^,Vr.), and since 

|^^+'-IdU,L,.+,<-Q-^<-l, (22) 
using Corollarv IA.2I we have the estimate 
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Now {gj,lvi} = for 1 < i < j from the induction hypothesis, and obviously we have 
{[f]vi,...,vjjvi} = for 1 < i < J, hence {ujjy.} = for 1 < i < j. Then if we apply 
Proposition 14.41 with 1^ = Ivi foi' 1 < ^ < J, we obtain that {u'^^l^.} = for 1 < i < j. But 
WjiUj+i) ~ ^-iicl hence {u'p /^j} = for 1 < i < j + 1. Moreover, we claim that {gj-,lvi} = 
for I <i < j implies {[gj]vj+i,lvi] = for 1 < i < j, and as {[gj]vj+i,lvj+i] = 0, this gives 
{[gi\vj+i-,lvi} = for 1 < i < J + 1, and this eventually gives {gj+i, Ivi} = for 1 < -i < j + 1. 
To prove the claim, note that since {lvi-,lvj+i} = for 1 <i < j, = ° -^Vj+i ^^'^ 

1 

s^+oo s In 



hm - r{gjoXi^^^,l,^oXi^^Jdt 

-4- + CXD S Jq 



s 

= 0. 

The estimate for gj+i is obvious: 

\gj+l\a,L,+i < \[gj]v,+i\a,Lj + \u'j\a,Lj+i < eQ'^ . 

Finally, from Lemma lA.ll and ()22p , we have 

\fj o ^j+i\a,L,+i < \ fj\a,Lj 

and therefore 

|/i+lU,L,+i < \Uj\a,L,+i + \ fj\a,L, <-eexp(- • Q'^^"). 

This proves the claim, and ends the proof of the proposition. □ 

4. Let us now prove Theorem II .31, which is an easy consequence of Theorem ll.il As before, 
this will be deduced from the more flexible proposition below. 



Proposition 4.6. Let H be as in ()G2p . and Q > 1. If 

Q->1, e<r^, r<A^{Q)-\ r < R, (23) 
then there exists a symplectic map $ G G°^{'Dr/2,'C'r), such that 

Ho<^ = h + [f]^ + g + f, {g,L} = 0. 
Moreover, for any I = (/i,/2) G N^", we have the estimates 

and 

\9L,L,i <• rVH'^lQ-\ \f\^^, <• rV-l'^l exp (- • Q-V") . 
Proof of Theorem \1.3l We choose 



with a well-chosen implicit constant so that the third part of (j23p is satisfied. Proposition 14.61 
with this value of Q implies Theorem 11.31 as the other conditions of are satisfied by dH) . 

□ 
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Now let us prove Proposition 14.61 

Proof of Proposition \4-6[ To analyze our Hamiltonian H in the domain Dj. , which is a neigh- 
bourhood of size r around the origin, we rescale the action variables using the map 

a:{9,l)^{e,rl) 

which sends the domain Di onto Dr, the latter being included in Dji by the last part of (I23p . 
Let 

H' = r-^{Hoa) 
be the rescaled Hamiltonian, so H' is defined on Di and reads 

H'{e, I) = r-^H{e, rl) = r'^h{rl) + r"V(6', rl), {6, 1) £Vi. 

Without loss of generality, we assume that h{0) = 0. Now using Taylor's formula we can 
expand h around the origin to obtain 

h{rl) =ru}-I + r'^ [ (1 - t)V^h{trI){I , I)dt = ruj ■ I + r'^h'{I) 
Jo 

and so we can write 

H' = L + f 

with 

f = rh' + r'\foa). 

Now we know that \f\a,L < £, so that the |r~^(/ o a)\a,L < r~^e, and using the second 
part of (I23D, \r-^{f o fT)!^, /^ ^ r. Moreover, applying Lemma IA.3I (with s — 2) we have 
\h'\a,L/2 <■ \h\a,L < 1, SO that \rh'\^^i/2 <■''• This eventually gives \ f'\a,L/2 <"^- Now we will 
apply Proposition 14.51 to the Hamiltonian H' = + f , defined on the domain Vi and such 
that |/'U,L/2 <■ ^- The first and third part of condition ()23p imply condition (|2ip with 1 
instead of R, L/2 instead of L and r instead of e, so that there exists a symplectic map 
G G"'^(2?i/2,Pi), with L = -L, such that 

H' o<^' = l^ + [f% + g' + f\ {5',U = 

with the estimates 

W-ld\^j^<Q-\ \g'\^;L<rQ-\ |/'^^ <• r exp (- • Q-V") . 

Note that [f']uj = [rh' + r^^{f o a)]^^ = rh' + r^'^[f]^ o a so the transformed Hamiltonian can 
be written as 

H'o<^' = l^ + rh' + r-'[f]^ oa + g' + /'. 

Now scaling back to our original coordinates, we define ^ = ao^' oa~^ , therefore ^ : — ^ 
Vr and 

Ho^ = rH'o^'oa~^ 

= r{L+rh' + r-'[f]^oa + g' + f')oa-' 

= {rL +r\')oa-^ + [f]^+rg' oa-^ + rf' o a~\ 
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Observe that (rl^ + r^/i') o a ^ = /i, so we may set 

g = rg oa , j =rf oa , 

and write 

Ho'^ = h+ [f]^ + g + f. 

The equahty {g,luj} = is obvious. Now let I = (Zi,/2) £ N^". Observe that from the 
definition of a and g, 

d^g = dl^d\^g = rr-\^^\d^g') o a'^ 

so 

|5'5lcop./,)<rr-l'^l|9VlcO(p,/,) 

and therefore 

Replacing (7' by /' in the above argument, we obtain 

l/L,z,/ < ^^""''l/'L,z,^ <• (- • Q"'^") • 

Finally, writing <I> = {(^0,(^1) and = ($g,<l>'^), we observe that (^e = ^'0 ° o""^ and = 
r^'joa~^ which immediately gives the estimates for ^0 and <!>/. This concludes the proof. □ 

A Technical estimates 

In this section we recall some technical estimates concerning Gevrey functions, taken from 
|MS02] . that are used in the proofs. 

Note that our Gevrey norm differs from the one in jMS02j . where they used 

instead of However these norms are "almost equivalent" in the sense that obviously we 

have \H\Qa,L(^jyR) ^ ll-f^llG^-^COfl) while, for instance, ||-ff||ga.L/2(-pH) — 2°(2"~l)~"^l-f^lG"-^(X'H)- 
So this change only affects the implicit constants. 

First recall that our proof of Theorem 11.11 proceed by induction on 1 < j < d, starting 
with the case d = 1 which is exactly Proposition 14. 3[ and the transformation in Theorem ll.il 
is a composition of transformations given by Proposition 14.31 Therefore we shall have to 
estimate the Gevrey norm of a composition of functions, and unlike classical norms used 
in the analytic case, these estimates are not completely trivial. 

Lemma A.l. There exists a constant C depending on n, a and L such that for any < 5 < R, 
if F e and ^> G C^'^' {Vr, ,Vr), with L' = CL and R' = R - 5, and if 

\^ - Idla^i/ <• I, then F o $ G G"^'^' {V ri and we have the estimate 
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This lemma is contained in the statement of Corohary A.l, Appendix A. 2, in |MS02j . 
in the case where is a real-valued function, but it extends immediately to vector-valued 
functions, by applying it components by components. The implicit constant also depends on 
d, but the statement will be used for a value of 6 depending only on d and R. 

Here's a direct consequence that we will also use. 

Corollary A. 2. There exists a constant C depending on n,a and L such that for any < 
S <R, if^ e and $ G C"'^' {Vr, ,Vr), with L' = CL and R' = R-6, and if 

1^ - Idla.L' <• 1; then ^ o <i) G G"'^'(P_r', M^") and we have the estimate 

1^ o $ - ldU,L' < 1^ - IdU,L + 1$ - IdU,L.. 

Let us also state a straightforward lemma which says that the Gevrey norm of a function 
controls the Gevrey norm of its derivatives, provided we restrict the parameter L (in the 
statement below, we simply choose L/l). 

Lemma A. 3. Let f G G'^'^{Vr). For any s G N, there exists a constant c depending on s,a 
and L such that 

sup |9'/U,L/2 < c\f\a,L- 

;gn2", \i\=s 
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